Neutral weak currents in nucleon superfluid Fermi liquids: 
Larkin-Migdal and Leggett approaches 
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Fermi liquids at finite temperatures. We explain peculiarities of both approaches and explicitly 
demonstrate that they lead to the same expression for the emissivity in pair breaking and formation 
processes. 
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Cn ' I. INTRODUCTION 

One of important mechanisms for cooling of superfluid neutron star interiors is nucleon Cooper pair breaking and 
formation (PBF) with a radiation of neutrino-antineutrino pairs, " N" — > " N" + v + v (see [ll— [20L [221 ] and references 
therein). Neutrinos are produced in weak interactions in which the lepton current — v(\ — 75)^ is coupled to 
the vector (V^ = <?y#7 /1 ^) and the axial-vector (A^ = <?a ^ 7^75^) nucleon currents; C — — ^) 

where G w 1,2 x 10~ 5 GeV -2 is the weak interaction coupling, gy = —1 and gA — 1.26. For the nuclear vector 
, ^ , ■ current holds d^V^ — that corresponds to the conservation of the baryon charge. Early works fll-fl4| did not 
care about the vector-current conservation. The latter is fulfilled only if in-medium renormalization of the vector 
7-H ■ current isperformed together with a corresponding renormalization of Green's functions. This problem was tackled 
in Refs. (T5rl2l| . Reference [l6| indicated that the emissivity of the ISq PBF processes on the vector current should 
be dramatically suppressed (oc vp, where v-p is the Fermi velocity of non-relativistic nucleons) provided the vector 
current conservation is fulfilled. Reference [1 61 ) used expressions derived within the standard BCS formalism of the 
superconductivity theory [23l. [2~H for low excitation energies ui, |g] <C A, where A is the nucleon pairing gap and uj is 
the net vv energy, whereas the PBF reaction kinematics permits only w > 2A, |cf| < oj. The correlation effects in the 
particle-hole channel were neglected and processes induced by the axial-vector current were disregarded. 

The consistent calculation of the PBF emissivity induced by the vector and axial-vector currents was performed 
in Ref. [ijj] within the Larkin-Migdal-Leggett Fermi liquid approach. The latter takes properly into account of 
the correlation effects in both particle-particle and particle- hole channels. It was demonstrated that the neutrino 
emissivity is actually controlled by the axial- vector current and is suppressed only by the factor tx vp, rather than 
oc vp as was stated in Ref. [Tfj| . This result was supported in the subsequent work which however continued to 
work out the vector current contribution neglecting the correlation effects in the particle-hole channel. 

The convenient Nambu-Gorkov formalism developed for the description of metallic superconductors, cf. Refs. [23l . 
[2^ |. does not distinguish interactions in particle- particle and particle- hole channels. These interactions can be, how- 
ever, essentially different in strongly interacting matter, like in nuclear matter and in liquid He 3 . The adequate 
methods for Fermi liquids with pairing were developed for zero temperature by Larkin and Migdal in Ref. [25] (see 
also [1^1) and for a finite temperature by Leggett in Ref. [13, [13 ■ The problem of calculation of a response function of 
a Fermi system to an external interaction becomes tractable at cost of introduction of a set of Landau- Migdal param- 
eters for quasiparticle interactions. Parameters can be either evaluated microscopically or extracted from analysis of 
experimental data, see [26]. The technical difference of the mentioned approaches is that Larkin and Migdal worked 
out equations for full in-medium vertices, whereas Leggett calculated directly a response function. The former ap- 
proach was aimed at the study of transitions in nuclei, and the latter on the analyzes of collective modes in superfluid 
Fermi liquid. The principal equivalence of both approaches was emphasized by Leggett in Ref. [13, HH . 

In Ref. [Uj] we used the Larkin-Migdal approach. More specifically, we solved the Larkin-Migdal equations for the 
vertices induced by the weak vector and axial-vector currents and calculated the current-current correlation function 
Xa, a = and the PBF emissivity. The explicit expression for Xa was obtained at T = 0. It is sufficient to 

calculate emissivity of the process for T <C 2A since small exponential factor e ~ 2A / T comes already from the phase 
space volume and the temperature correction to Xa(T = 0) is small in this limit. Despite Larkin-Migdal equations 
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FIG. 1: Amplitudes of the transition of a particle into a hole and a condensate pair, A^ 1 ' and vice versa, A^ 2 ', and normal, 
G, and anomalous, F, Green's functions. 



were derived in their original paper for T = 0, actually the results can be generalized for arbitrary T. In Ref. [18j we 
sketched which expressions should be modified at finite temperatures. 

Recent paper [22j tried to adopt the Leggett formalism to calculate the PBF emissivity for the ISo neutron pairing 
and came to different results compared to those derived in Ref. [Hj], even in case T <C A. One of the points of 
Ref. [22[ was to find the PBF emissivity for arbitrary T ^ 0. In view of the explicit claim by Leggett in Ref. [28[ 
on the equivalence of his approach to that of Larkin and Migdal, this difference looks worrisome and requires a 
clarification. 

The aim of this paper is to reveal the correspondence between the Larkin-Migdal and Leggett approaches and to 
generalize the results of Ref. 18 1 to arbitrary temperatures. We argue also that the results of the work [22| are based 
on misinterpretation and wrong solution of the Leggett equations. 

The paper is organized as follows. In Section [TT] we introduce the Fermi liquid approach to the problem of the 
neutrino emissivity from nucleon matter. We focus on the emissivity via the PBF reactions. The main expressions 
are presented in the diagrammatic form valid at arbitrary temperatures. Then in Section [IIII we demonstrate how the 
Larkin-Migdal equations and the Leggett equations follow from the same set of diagrams. In Section HVl we solve the 
Larkin-Migdal equations for the vertices induced by the vector and axial- vector currents and, then, present expressions 
for the PBF emissivity recovering the results of Ref. 18[ formulated now for arbitrary temperatures. In Section IVl we 
demonstrate how the same results can be obtained within the Leggett formalism. In Section I VII we discuss the flaws 
in Ref. [11. We conclude with Section [VE] 



II. GENERAL EXPRESSIONS FOR NEUTRINO EMISSIVITY AND CURRENT-CURRENT 

CORRELATORS 

In this section we recall the main expressions for the neutrino emissivity obtained in Ref. [18], writing them in 
the form valid at arbitrary temperatures. To be specific we consider a neutron system with the s-wave pairing. The 
neutrino emissivity 

-— f TT d3qt ( a;i + ^ 2 ) s Sx(gi + m 

£ ™ 4 J 11 (27T) 3 2Wi e (a»i4**)/T _ ! W 

is expressed through the imaginary part of the Fourier-transform of the current-current correlator x(l) = Xvil) + 
XA ( q ) = J d 4 xe-^ q -^(N\(V^(x) V v (0) + A»(x) A v (0)) l^(x)lt(0)\N) , where q = (wi + u 2 ,qi The averaging is 

done over the vector of state, \N), of the fermion system with pairing at thermal equilibrium. The summation runs 
over the lepton spins. 

We use non-relativistic limit for nucleons since the nucleon Fermi energy €f *C m, where m is the nucleon mass. The 
bare vertices for the nucleon currents as they follow from the expansion of the Lagrangian are V 0tl rj ip^ Op') (1, (p + 
p)/2m) ip(p) and A ^ ps tp^(p') (<?(p' +p)/2m, a) ip{p) , where a = (01,02,03) are the Pauli matrices acting on free 
nucleon spinors ip, and p and p are outgoing and incoming momenta. 

In a system with pairing a particle can transit into a hole and a condensate pair and vice versa. The one-particle 
one-hole irreducible amplitudes of such processes are depicted in Fig. Q] Besides the normal Green's function G 
one introduces anomalous Green's functions , and F^ , which diagrammatic notations are given in Fig. [TJ The 
Green's function for the hole is defined as G h (p) = G T (—p) . Superscript T denotes matrix transposition. The normal 
and anomalous Green's functions are related by Gor'kov equations. The normal Green's function is diagonal in spin 
space G(p) = G(p) 1, but the anomalous Green's functions are equal to F^(p) = F^ jp) = F(p ) = F(p)i0%. We 
will not discuss here the details of the Fermi liquid renormalization procedure, see Refs. [25H27;. 32]. The latter allows 
to separate contributions from the pole and regular parts of the Green's functions in all relevant quantities. We 
assume that this procedure is properly done and we deal further only with the pole parts of the Green's functions 
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FIG. 2: Diagrams contributing to the susceptibility \ m ©• 

characterized by the effective mass to* and the residue a: 

C( \ = a ( £ + e p) F( \ = ~ aA (9 s 

[P> e 2 -E^ + iOsgne 1 {P> e 2 - E 2 + iOsgne ' l j 

where i? 2 = e 2 + A 2 with the pairing gap A, and e p — p 2 / (2 to*) — /i(T), /U is the chemical potential, /i(T = 0) = ep = 
Pp I (2 m*). At finite temperature one can use the Matsubara techniques with the replacement e — » % e n — i 7r(2n + I) . 
The amplitude in the particle-particle channel is parameterized as 

[ f? ]L C = T o(n,n')(ia 2 )t(ia 2 ) c d + rf(n, n')(ia 2 d)t{cr i<r 2 )% 

and the interaction in the particle-hole channel is 

[f "] a b c d = T%(n, n') 51 8 d + rftrf, n') (a) a b (a f d . 

Here and below n = | and n = p'/l^'l • Superscript "cj" indicates that the amplitude is taken for \qv f| <C w 
and w <C £f: where and q are transferred energy and momentum. We use that ~ e p + vq/2, where 

V is the nucleon velocity at the Fermi surface, v = vfH(1 + 0(T 2 /e F )). Since in the PBF processes \q\ < A, 
lu ~ 2A, the terms neglected are of the order of 0(A 2 /p F ). Actually, the denominators of the Green's functions are 
\u ± (fp+g/2 — cp*-g*/2)l ^ e F- Moreover, the terms oc vq may vanish under the angular integrations. Taking this into 
account we estimate that the neglected terms are at most of the order of A/ep <C 1 compared to the retained terms. 
Such corrections are usually omitted in most calculations within the Fermi liquid theory for superfluids. 

The empirical information is usually expressed in terms of dimensionless parameters 

where p — m ^i F is the density of states at the Fermi surface. 

The pairing gap is determined by the Tg term in the particle-particle interaction, and the gap equation reads 

A(n) = ft) MW)) A(n'))n' , (4) 

A (A) = J d<P T G ( P )G h (p)6(Z-e p ), (5) 

where Gq(j>) = l/(e — e p + iOsgne) is the Green's function for the Fermi system without pairing (A = 0), and £ ~ /i 
is the cutting parameter. Usually the gap is determined with £ = fi. We use the following notations for the angular 
integration 

^(..•), (6) 



and for the integration 



r°° de r°° 

>L W-„ de ' /fce » )forT=0 

pT ^2 / de )3 /(ie n ,ep) for T ^ 
Note that the same value A can be introduced as A = f d$ T (G(p) G h (p) + F(p) F(p)) 0(£ - e p ), cf. Ref. [11]. 1 



/" d$ T /(e,e P 



1 Definition of the value Aq is here the same as in Ref. |25| and differs by sign from that used in Ref. [27, 28]. 
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FIG. 3: Graphical representation of dressed vertices in Fig. [2] 



The diagrammatic representation for the current-current correlator entering ([]} is shown in Fig. [2] The dash line 
relates to the Z-boson coupled to the neutral lepton currents. In terms of the Green's functions and the quasiparticle 
interactions we write it as 

Xa (uj, q) = ^Tr( J d<t> T rZ { G+ f a G_ + F + r h a F_ + G+ f«F_ + F + t< 2 > G_ } ) , (7) 

a = V for the vector current and a = A for the axial-vector currents. The trace is taken in the nucleon spin space. 
Doing calculations in Matsubara technique, we use notations: G+ = G(po + uj,p + q/2) and G— = G(po 1 p — <f/2), and 
similarly for F± functions. For continues frequencies we use the symmetrical 4-vector notations, i.e. G± — G(p±q/2) 
and, analogously, F± — F(p ±q/2). The left vertices in Fig. [2] are the bare vertices generated by the weak currents 

Ty(n, q) = 9V (Tv,0 lo-Tv,l_l)i, Ty,0 = , fy A = ~~"' ( 8a ) 

t%{ft,q) = -g A (TA,ivlo-TAflV l )> T lo = ^-, f '^' 1 " : ^" 17 ' ( 8b ) 

ey and e A are effective charges of the vector and axial- vector currents, cf. [18j . The corresponding vertices for holes 
are defined as 

Tv' h (n,q) = [T^(-ri,q)] T = g v (Tv, lo + Ty A l)l, (9a) 
(fl,q) = [TA(-n,q)] T = -9A ( ~ UK7 T i - tXo<t T Z) . (9b) 

We used here explicitly that T% (-p,q) = 7% i0 (p,q) and f"^ A (-p,q) = -T%,i(p,q)- 

The right vertices in Fig. @ are the full in-medium-dressed vertices, which are functions of the out-going frequency 
w and momentum q, and the nucleon velocity v ~ vf n, ft — p/p . In absence of external magnetic field the spin 
structure of the full vertex is the same as the spin structure of the bare vertex. Therefore, we write 



tv = 9v (rv,o lo ~ tv,i l) i , Ty = g v (t v>0 1 - f y,i I ) 1 , (10a) 

4 1} = (4 2) ) f =-9v(rv,olo-T V ,J)ia 2 , (10b) 

ta = -g A (ta,i&Io - TAflvl) , t\ — ~9a (t a,i^ t h - tajq ) , (10c) 

T A = ( t a)^ = +9 A {j A ,ldlQ - TA,0&1) i<?2 ■ (10d) 



We anticipate here the relation between the vertices f' 1 ' and t^ 2 \ which will be proven later. The full vertices are 
determined by the diagramatic equations depicted in Fig. [3] These graphical equations were first introduced by Larkin 
and Migdal [25j . The blocks in Fig. [3] correspond to the two-particle irreducible interaction in the particle-particle 
channel, T^, and the particle- hole irreducible interaction in the particle- hole channel, T u . We emphasize that only 
chains of bubble diagrams are summed up in this particular formulation. Thus, the imaginary part of x accounts 
only for one-nucleon processes. To include two-nucleon processes within a quasi-particle approximation, one should 
add diagrams with self-energy insertions to the Green's functions and iterate the Landau-Migdal amplitudes T^'* in 
Fig. |3] in the horizontal channel (29j . 

Working out the spin structure in (JT]) we obtain 

Xv(q) = 9v(( T v,o l o-T V il) (l xv,o(n,q) - Xv,i(n,q)l ))-, (11a) 
XA{q) =9a((ta,iIo -Ta, 1) {ioXA,i(n,q) - XA,o(n,q)l )) n , (lib) 
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where we introduced new scalar and vector response functions 

Xa,o(n, q) = J d$ T (G+ G_ r Q , - F + F_ r£ + (G+ F_ - F + G_) r O)0 ) , (12a) 

&,ifa«)= [d*T(G + G-?aA+F+F-T* 1 + (G + F--F + G-)? a ,o)- (12b) 
It is instructive to express susceptibilities (jllaj) and (|llb() in terms of the polarization tensors 

x (?) = W(<?)^. 

The polarization tensors can be written through the auxiliary bare and dressed currents, Jau{^Q) anc ^ Ja^ffi, q), 
respectively, 

K V (<1) = \(Ti-{J^(n,q) r a (n,q)})n , (13) 

where the vector currents are 

Jv ^(n,q) = (iv,o,Tv,i) 1 , Jyi^l) = (xv,o?xV,i) ^ ( 14 ) 
and the axial- vector currents are 

J^i^l) = ( f?f *A,i' (?r Xo) > Jli^Q) = (fXA^VXAfi) ■ (15) 
Thus, the polarization tensors can be cast as 

^ _ ( D$' A _ / (7t?, Xv,o)s, W,o(fy,i)i)s , rifrt 



The integration over the lepton phase space can be performed analytically, cf. [181 ] . and the neutrino emissivity is 
then cast in terms of the polarization tensor as 

C 2 f°° f u tin 2 
' 1 dc / d \q\-J—K a (q), (18) 



48 ^ Jo Jo |y 

K a (q) = {q^q v -g^q 2 )^{q). (19) 

The polarization tensor for the conserving vector current must be transverse Hy^q^ = Tly' = 0. This property will 
be explicitly proven below in Section IIVI and in Appendix |B] Taking it into account and using (|16|17j) we find 

Ky(q) = ^A(q 2 - L o 2 )^( Xv , (n,q)~-vxv,l(n,q))n, (20) 
a 

a 

- u(qx AA (n,q)) ri -oj((qv)xA,o( ri >q))fi\- (21) 

Finally, once the full in-medium vertices (|10all0bll0cll0dl) are known, expressions (I18I20I21[) together with (|12all2b|) 
solve the problem of neutrino emission from superfluid neutron matter via the PBF reactions. 
In the rest of the paper we demonstrate two methods for the calculation of the full vertices. 



K A{q) = -^-^3[g 2 (w£-u(^,<7))n + (3w 2 - 2q 2 ) (xAfi(n,q))n 



III. CORRESPONDENCE BETWEEN LARKIN-MIGDAL AND LEGGETT FORMALISMS 

The coupling of an external field to the non-relativistic fermion is described by the 2x2 matrix acting in the fermion 
spin space. Any rank- 2 matrix can be decomposed into a unity matrix and Pauli matrices a. Thus, we have 

f(n,q) = t (n,q)l + ati(n,q) , 

f«(n,g) - (t 1 \n,q)l + at 1 il) (n,q))i ( r 2 , 

f ^ (ft, q) = ia 2 (4 2) (ft, q)l + a f« (ft, q)) . (22) 
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The hole-vertex is decomposed as 



(ft,q) = tZ(n,q) 1 + a t?(n, q) = t (-n,q) l + aH(-n,q) 



(23) 



The spin structures of the weak coupling vertices (O and (|10[) , we primarily deal with, demonstrate that the vector 
current contributes only to the vertices with subscript and the axial- vector current couples only to the vertices with 
subscript 1. E.g., for the vector current 



and for the axial-vector current 



'o 


= 9v 


\ T V.O k 


to 


= 9v 


[t~v,o k 


ti 


= h = 


= 0, 


to 


= to = 


o, 


tl 


= -9A 


\ T A,1 


tl 


= -9A 


(ta,i 



(24) 



- ta,o l) ■ (25) 

After opening the spin structure of the diagrams in Fig. [3] we arrive at the following set of equations for to, ig, t^p 
and t (for brevity we omit the dependence of the vertices on n, oj and q ): 



to — 1 — 



d$ T G+ G_ t - F + F_ <S - G+ F_ t 1} - F+ G_ # 



,(2) 



l C 



d$ T rg 



G h + G h _ t' Q l - F+ F_ t - F+ G h _ 4 1} - G\ F_ i 2) 



(1) _ 
o — 

(2) _ 



J d$ T T« [g+ G h _ 4 1} - F_ 4 2) + G+ F_ t + F+ G h _ )^ 
d$ T T« [(& G_ 4 2) - F + f - *o 1J + F + G - to + G+ F- tl 



(26a) 
(26b) 
(26c) 
(26d) 



The similar set of equations for 3- vector vertices t\, ti, and t^' is written with the only differences that Tq^ 
is replaced by and in front of all terms with t h the sign must be changed. The origin of this sign change is the 
identity 02 <7 T 02 = —a . 

For the sake of convenience we introduce brief notations, e.g., 



d$ T G h , F_ =Gi-F_. 



(27) 



The details of calculations of these products within the Matsubara techniques are deferred to Appendix [XJ E.g., we 
recover useful relations 



G f~ih f~ih /~i si Tp 771 fi zp f~ih s~ih rp 

-4-'Lr_ — Li ' I i , Lt_[_-_T_ — — _T_|_-Lt_ , — — Lr _j_ ■ _T _ , 



(28) 



see Eqs. (|A13|) . (|A15|) and (|A18[) in 28]. From (|2"B1) we can immediately find relations between the vertices tg and 
t 2) . Taking the sum of Eqs. (|26c[) and (|26d[) and making use of Eq. (|2"51) we obtain the homogeneous equation for the 
sum 4 15 +*o 2) ) 



L ' L 



d$ T r« 



G+ Gl - F+ F_ 



which implies 



(l) +t (2) _ Q 
f -|- I — U . 



(29) 



(30) 



The latter relation justifies the parameterization of the full in- medium vertices in Eqs. (|10b[) and (|10d[) . The same 
relation is valid for t \ ^ and t \ ^ vertices. 
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A. Larkin-Migdal equations 

In their original paper [25[ Larkin and Migdal presented Eqs. (l26l) in somewhat different form. They noted that 
the vertices for the holes, t h , can be obtained from the particle vertices with the change n — > —ri, 

to (n,q)= t (-n, q) , t 1 (n,q) = t 1 (-fi,q). (31) 
Therefore, one can introduce the operator tt, which performs this change of ft in the vertex 

TTt (n,q) = to(n,q) , nt x {n,q) = t\{n,q) . (32) 
In view of relation (|30|). Eqs. (|26cl26dl) reduce to one equation for the vertex 



Analogously we introduce t\ — —t^p — t\ ■ Then four Eqs. (|26p for scalar vertices '0' and four equations for 3-vector 
vertices t\ can be cast in terms of four equations 



to -US = (r£f(£(7r)to + M*,))„, (33a) 

t = -(r^({N + Ao)ta + 0(it)to)) , (33b) 

ti-tj = (ry^-^fi + Mfi)) , (33c) 

ti = -(rf (^(N + A )t 1 +O(-7r)t 1 )^ / . (33d) 



We shall call this set of equations the Larkin-Migdal equations. In Ref. [25] these four equations are further reduced 
to only two equations with the help of the operator P (defined by Eq. (31) in Ref. [25|). which includes additionally 
the change of the sign between Eqs. (|33aj) and (|33cj) and between Eqs. (|33bj) and (|33dj) . Functions L, M, N, and O 
are defined as in Ref. [25[ 

L(H,q;n) = G+ ■ G_ - F+ ■ F_ tt , 

M(n,q) = G+-F_ - F+-G- , 

N(n,q) = G + -G h _+F+-F--A , 
0(n,g;7r) = -G+- F_ - F+-G h _K ■ (34) 

We emphasize that Eqs. (l33l) are valid at arbitrary temperature. The temperature dependence is hidden in the 
convolutions of the Green's functions (|3"4"|) . In Ref. (25| the latter ones were calculated explicitly only for T = using 
the method of Ref. 30]. The extension to T ^ is straightforward within the Matsubara technique, see Appendix [Al 

B. Leggett equations for vertices and correlation functions 

In this section, starting from Eqs. (I2"r?|) . we derive the equations obtained by Leggett in Refs. [13, [28|, thus, demon- 
strating interrelation of Larkin-Mi gda l and Leggett approaches and their principal equivalence. 

The primary goal of works (27l |28{ was to calculate the response function of the superfiuid Fermi liquid to the 
source of an external field f"(n, q) at uj, q <C A and to study collective modes in a superfiuid. The excitations of the 
density, spin-density, current and spin current fluctuations with f w = 1, 03, n and ft 173, respectively, were considered. 
In Refs. [27], [H[ it was explicitly assumed that the source is diagonal in the spin space and the vertex is characterized 
by the directions of n and of the projection of the spin s. Then the vertices can be classified as "even" and "odd" 
depending on whether they change the sign at the simultaneous replacement ft — > —n and s — > — s . So the vertices 1 
and n (73 are " even" , whereas 173 and n are " odd" . 
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Using relations (|28| we are able to present Eqs. (|26l) in the matrix form 



%(n,q) = t^(n,q) + (f (n,fi') G±(n' , q) %{n', q)) n > 



'■q 

f (2) 



'0 

to 

-G% ■ G- 
-F+ ■ F_ 
G\ ■ F_ 
F+ ■ G- 





■ " 













, r = 


' '0 _ 


j-hoj 
l 

- z - 











o r« o o 



o o 
o o 



o n 



-F+ ■ F- G\- F_ F+ ■ G- 
-G+-G h _ -F+-G h _ -G+-F_ 
-F + -G h _ G\ ■ G h _ -F+ ■ F_ 



-G+ • F_ 



-F + ■ F_ 



G~\- ■ G— 



(35) 



(36) 



The matrix Ga is precisely the matrix used by Leggett in Eq. (3) of Ref. [281 ] . The same matrix enters the equation 
for the vector vertices in 



fi(H,q) = f"(ri,q) + {Tx(ft,fi')Gi(n , ,q)fx(n , ,q))r{i , 



(37) 



Tx = 



r r (1) " 
1 1 




' 




1 1 













, fx = 


-h 
-tx 


> >1 ~ 






. tx . 




-KjJ 

t 





r; 

rf 

v% 





These equations elucidate the meaning of the "arrow space" introduced by Leggett in Ref. [27[ in application to the 
in-medium vertices. 

To proceed further we introduce even and odd vertices. In our case the even vertices are 



and the odd ones are 



t + — + *o) > t + — — t{ 1 ) 
t = |(to — t ) > t = |(ti + tx) 



(38) 
(39) 



With these definitions the bare vertices 1 and na^, considered by Leggett in Ref. [28|, contribute only to the even 
vertices, whereas ft and 03 to the odd vertices. 

Equations for the even and odd amplitudes follow from Eqs. (|26l) if we take the half-sum and the half-difference of 
Eqs. (|26a|26b|) . Taking the half-difference of Eqs. (|26cl26dj) and making use of Eq. (|30[) we obtain the equation for 
the anomalous vertex ^(t^ — t^) — — t. The systems (|35|) and (|37|) acquire then the same matrix form 



t(n,q) = T"{n,q) + (f (n, n)G 3 (n' , q) T (n ', q)) H , 



(40) 





' -t 









"r* " 


t = 


-t~ 


, f" = 


-ir u 


, f = 


r w 




t+ 








r w 



-{G+-G h _ + F+-F_) G h + -F_-F + -G- 
G\-F- -F+-G- \{G+-G- + G h + -G h _) + F+-F- 
G h + -F_+F+-G- \(G h + -G h _-G+-G-) 



G\-F- +F+-G- 
\(G h + -G h _-G+-G- 
\(G+-G 



G h + -G h _) - F+-F_ 



(41) 



Lq (?) for the 



and we do not distinguish here the vertices with subscripts and 1. In Eq. (j40|) we have to use T u< & 
vertices without a matrices (subscript 0) and L u< ^) = for the vertices with a matrices (subscript 1). In Eq. fH 

we recognize the minor of the matrix g introduced by Leggett, see Eq. (12) in Ref. |28l|. Actually, Leggett in Ref. [25 
presented these equations in a different form. To reproduce this form, following [28J we introduce new quantities 



K (n,q) = \{G h + -G h _ + G+-G-) + F+ • F_ , X(n,q) = F+-F_ 



(42) 
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All entries in the matrix (l4"TT) can be expressed through these two functions 

G% ■ F- - F+ ■ G_ 
G h + -F-+F + - G_ 

G+-G h _+F+- F_ 
\{G h + -G h _-G + -G_) 



^A(r?,g) , 



= An 



{vqf 



2 A 2 



X(n,q), 



-—K{n,q). 
v q 



(43a) 
(43b) 

(43c) 
(43d) 



In |28| Leggett emphasized that relations (^a^b^cl) ) are valid for any uj, q and T, whereas relation (|43d|) holds 
only in the limit vf \q\, to -C A. In Appendix lAl we re-derive these relations and demonstrate that, actually, Eq. (|43d|) 
is valid for arbitrary u) and q. Thereby, we prove that one may use these conditions in kinematic region of the PBF 
reactions. 

With the help of Eq. (|4"3"]l we rewrite the system of equations (|4Tfl) as 



" -t 









"r* o o " 




-f- 




-t- u 




o r w o 




t+ 








o o r u 





-Ac 



" -{v<D \ _1L2 \ 

_2A 2 A A 
vq y 

~~A 
A A 



A 



A 



- — K 

v q 

k-2X 

v q 





' -t 




-t~ 




t+ 



(44) 



-t -> V>2 , * + -> V>3 ■ 



(45) 



Here we recognize Eq. (22) of Ref. [28[, if we identify 

t +u -> e + , -> c , -* -> ^i 

Note that according to these assignments the quantities ^ ± introduced in (22) by Leggett ought to be defined as 
^(n) = |[£(— p, —s) ± £(p, s)] . Leggett defined £~ with opposite sign. It has no influence on the correctness of his 
results, since the quantities V'1,2,3 were not identified in |28j with in-medium vertices and the final expressions for 
the correlation functions, Eq. (23) in Ref. [28[ depend quadratically on ^ . Further, we shall name the system of 
equations (|44l) . the Leggett equations. 

Now we formulate the current-current correlation function in terms of the Leggett notations (using t and t vertices). 
First, we express correlation functions ([7]) in terms of vertices t, t h and t, cf. Eq. (|2"2"|) . 

X = Xo + Xi . 

Xo = ( j d$ T US (n, q) (g+ G- t (n, <?) - F+ F_ tg (ft, q) + (G+ F_ -F+G-) to (n, q) 
Xi = ( / d$ T t?(ft, q) (G+ G_ h(n, q) + F+ F_ t> l (n, q) + (G+ F_ - F+ G_)? x (n, q)) ) r . 

Further transformations we illustrate at hand of \o- We do replacements e p — ¥ —e p and e n —¥ — e„ (e — > — e for 
T = 0), which do not change the integrals and the Matsubara sums, but induce replacements F± —> F T , G± O G^ 
and f O t' 1 . Then we add the resulting expression to the initial one, divide by two and obtain 

Xo = ^(Jd^ T t^(n,q)(G + G_t (n,q)-F + F_^(n,q) + (G + F_-F + G^t (n,q)) 
+ ±(J d^ T t^(n,q)(GlGt^(n,q)-F + F_t (n,q) + (F + G h _-GlF_)t (n,q) 
The latter expression can be easily rewritten in terms of the even and odd vertices tg with the result 
Xo = (4" [{G h + -F. + F + ■ G_ ) (-t ) + \ {G\- G h _ - G+ ■ GJ) (-<„ ) + ( | (G+ • G_ + G%- G h _ )-F + -F_)t+ 

+ ({-t^) [{G h + -F_ - F+-G_) (-t ) + (i(G+-G_ + G%-G h _) + F+-F_) (-t„ ) + \{G h + -G h _ - G+-G_)t+] )_ 
which in the Leggett 's notations (|4"2")) turns into 



Xo 


= Xo+ + 


Xo 




= (tf 




Xo+ 


A 


Xo- 


= ((-to 





w q 



v q 
A" 



A(-t ) + «(-t ) 



■ K tr 



v q 



(46) 
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With replacements (|45|) we recover Eqs. (23a), and (23b) of Ref. [28| , Note that the current-current correlator \ i s 
given by the sum of the even and odd terms, whereas Ref. [28j presents two independent expressions, one oc £ + and 
another one oc £~ . For the vertices considered in Ref. [28[ either the term oc £ + is zero or that oc £~ . But in general 
case the correct expression is given by Eq. (|46[) . Similar equation holds also for the vertices with subscript 1: 

Xi = Xi+ +X1-, 
Xi+ = *i - x X(-t 1 )-—K(-t 1 ) + ( K -2X)t 1 ) 

Xi- = Hi "x A Hi ) + «("*! >_■ 47 

\ LA v q 1 / fi 

Concluding this section we stress that both Larkin-Migdal and Leggett equations for the vertices and the current- 
current correlators follow from the same set of equations, and, hence, the results of the calculation of the emissivity in 
both approaches should be the same, provided calculations are performed correctly. Now we will focus on the solution 
of these equations. 

IV. SOLUTION OF LARKIN-MIGDAL EQUATIONS AND CORRELATION FUNCTIONS 

We apply the Larkin-Migdal equations (|33[) for the case of the weak-current vertices (|9ll0p . For the weak vector 
current vertices we use Eqs. (|33al33b[) and for the weak axial- vector current vertices, Eqs. (I33cl33d|) . Then we separate 
the parts proportional to the scalar Iq and to the vector I and obtain altogether 8 equations for vector and axial- vector 
current vertices. We cast these sets of equations in the following form [18| . 

T a ,o(H,q)=T^ (n,q) + (T^(n, n') [L(n<, q; V a . ) T a , (n< , q) + M(n', q) T a , (n', q)] ) R , , (48a) 

r„,o(n,<z) = - (ri(n,n') [{N[n>, q) + A ) ?„, (n', q) + 0(n>, q; V a ,o) r„ j0 (n', q)] , (48b) 

f aA (n,q) = < 1 (n,q) + (T^(n, n') [L(n> , q; V a ,l) T a ,i (n> , q) + M(n>, q) T a ,i (»', ?)] ) H , , (48c) 

r a>1 (n,g)= - (r«(n, n>) [{N{n>, q) + A ) ? a<1 (n', q) + 0(n', q; T a ,i) r a ,i(< ?)] (48d) 

To write one set of equations for both vector and axial-vector weak currents we introduced the notation for the 
effective interaction = Tq , if a = V, and T"^ = , if a = A. Operators V a ,i are defined as follows 

V a ,i = {-lfP a ^, i = 0,1 (49) 

with parameters 

P Vfl = l, JV,i = -l, P Afi = -l, Paa = 1, (50) 
which are eigenvalues of operators V a ,i, when the latter are acting on the bare vertices 

£a,o< = P a , r a V ^1=^,1^1- (5!) 

To proceed let us for simplicity assume that and T| contain only zero-th Legendre harmonics. Extension to 
higher harmonics will be done elsewhere [3lJ. From (|48b[) we find 

(0(n,q;P afi )) H 



r a ,o(q) = -ril /AT ,^ / x " r a ,o(q), (52) 



where 



e _ Ti{N(n,q)) H 



s _ - a r v_2 (53) 

' a l + rl^o + Wn,^]' 

For the channel a, for which the gap equation 1 = — Ao T| is valid, we obtain r/^ — 1. For another channel, rj^ ^ 1. 
For the s-wave pairing we consider in the present paper, it holds 1 = —A a Yy. In Ref. 18f\ we have put rj| = 1 in 
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both channels since owing to the identity (0(n',q; — 1))-, = one gets ta.q — independently on the assumed value 
of ri A , provided only zeroth harmonics of the interaction are retained. Substituting Eq. (|52p in Eq. (|48a[) we obtain 

7-0,0(9) =7a(g , ;io,o) t-o,q> 

1 - 1 (q;P) = l-T^(C a (n,q;P)) ii , (54) 

where we introduce the notation 

C a (ft, q; P) = L(n, q; P) - r£ ^ ?) " (55) 

Solving the second pair of the Larkin-Migdal equations (I48cl48d|) we first note that for the constant and rf the 
angular averages on the right-hand sides of equations do not depend on n. Therefore, the component of the bare 
vertex proportional to v is not renormalized in medium. However, in view of the identity 

(f(n,q)n) n = (f(n,q)(nq)) n q/q 2 (56) 

valid for an arbitrary scalar function / of n and q, the full vertices gain a component proportional to q. Thus we 
decompose 3-vectors r aj i(n, q) and r a ,i{n,q) into the parts proportional to the n, ri q = q/\q\ vectors and introduce 
new scalar form factors 

T a ,i{n,q) =r" A n + T^l(q)n q , T aA {n,q) =T ( a q l(q)n q (57) 

with the bare vertex = (nr " ± ) . Action of the operator V at % on the vertices (|57p is given by 

P a ,lT ai i(n,q) = P a ,iT^ 1 n+ {-P a ,i)n q T^l(q), 
Va,lT a ,l(n,q) = (-P atl )n q T { a %). 

Then, from Eq. ([48d|) we recover 



-(<?)_ £ (Q(n,q;-Pg,i))n (g) s (0(n,q;Pa,i)(nn q ))n : 



From Eq. (|48c|) . substituting there Eq. we find 



r { a q l = 7a(q;-P a ,i)K(C a (n,q-P aA )(n •n g ))^< 1 . (59) 

Here we introduce the quantity 

£ a (n, q; P) = L(rt, q; P) - 4^M^T °&> « F ) ( 60 ) 

and use the identity 

(£ a (n,q;P a ,i))n = (C a (n, q;P a ,i))n , (61) 

which allows to use for the vector vertices the same function j a as for the scalar vertex . 
In terms of the loop-functions (|34|) the response functions (|12all2b[) can be expressed as 

Xafi{n, q) = L(n, q; V a ,o) T afi (n, q) + M(n, q) r afi (n, q), 
Xa,i(n,q) = L(ri,q;Va,i)T a:1 (ri,q) + M(n,q) T a j.{n,q). 

Using solutions (I5"2"j) and ([54]) for the scalar vertices we find 

Xa,o{n,q) = j a {q;P afi )T^ C(n,q;P afi ) . (62) 

With the help of Eq. (1571) we construct 

XaAfi'l) = L (n,q; p a,i)fir"i +L(n,q;-P ai i)H q T^l(fi,q) + M(n,q) n q T ( a q l(n,q) . (63) 
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Using solutions (158)) and (|59j) for the three-vector vertices we obtain 

Xa,i{n,q) = [L(n,q;P a ,i)n- M(n,q)rjl — H - n q )T x 

+ C a {n,q]-P aA )^ a {q]~P a ^)T^{C a {n,q]P aA ){nn q )) ri n q T^ 1 , (64) 

and, then, rewrite it as follows 

Xa,i(n,q) = T^ :1 j a (q;-Pa,i) C a (n,q;P aA ) + 5xas(n,q) , (65) 

5Xa,i{n,q) = vt r^f^jy, (°( n '>^ P a.i){n - n'))-, t% a 

+ j a {q; —P a ,i) t£i [c a (n, q\ -P a ,i) (C a (n, q; P a ,i)n)n - C a (ft, q; P a ,i)n (£ a (n, q; —P a ,i))nJ ■ 

Please, pay attention to some misprints in Eqs. (16,17) of Ref. (l8| . which are now corrected in Eqs. (|58I59I62|) and 
(|6"5"1) . The final expressions for the emissivity in [18[ remain unchanged. 

In [l8T | we have demonstrated the vector current conservation up to terms oc Vp for T = 0. In the present work we 
are able to prove that the vector current is exactly conserved for arbitrary temperatures. In Appendix [B] we prove 
relations 

(^xv,o -qxvs)n = o, 

^((qv) (uxv,o - qxvs))n = 0. (66) 
These relations ensure the transversality of the polarization tensor for the weak vector current (|16[) . H^q„ = = 

ntiu 
... V ■ 

The averages entering the quantities Ky and Ka, Eqs. (|20p and (I2ip . which determine the neutrino emissivities 
(1181) , acquire now the following form 

(Xa,o{n,q))n = -fa(q;Pa,o)Ta, Q (£(n,q;Pa,o))n, (67a) 

(Xa,o(n,5) {qv))fi = j a (q;P afi )T^ Q (£(n,q;P a ,o)(qv))a, (67b) 
(qXa,i(n,q)) H = j a (q;-P aA )(£ a (n,q;P at i)(T"iq))n- (67c) 



(vXa,i(n,q))n = (L(H,q;P ail )(v r " >1 )) s - (M(n,q) (vH q ))nvi 

+ Ja(q; -P ,i) r a (£a(n,q; -Pas) (vft q ))fi (C a (n,q; Pa,i){Ta,in q ))n ■ (67d) 



(N(n,q))rt 



Let us now find the neutrino emissivity at the condition tip < 1, being valid in the region of the IS neutron pairing 
in neutron stars fJ-TTstl - First, we expand the averages of the £ and £ functions 

— (£ a (n,g;+l))« « — f- I 1 + J g T (0, 0, 0) + — — [ — + fl ,^ 2 



a p ouj \ bui 1 5w V ow 2 o[vq) 

^■^(O.w.Qj + ^-l) 

..2 a -<4 4 



_ flr (o lW> o) + (^-i)— ww - 



1,(1^)^5^0)^^^, ( 68 a) 



^22, o ->2 2 



( 1+ 5^) (5t(0,0,0)- 5t (0, W ,0)) 



3 a; 2 V 5aj 

4 



4<? 4 /3 5T (0,0,0) , 9gr(0,0,0) 



<9w 2 d(vq) 2 
g 4 4 dgr(0,^,0) 6 6 6 

5^ <W + ° (<Z ^ } ' (68b) 
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-5— {C a {n, q; +1) (nn q ))ft 
a A p 



{£a(n,q;-l) (nn q )) f 

ar p 



\q\v F 



4\q\ 3 4 

3lo V ' 5u) 2 4 5w 3 



1 



5t(0,^,0) 



v F q A fdg T (0,0,0) a ffT (o,o,o) 



5w 



d(vq) 2 



(^-U^fi + ^Wo^o) 



3w 

|3 ,,3 



3w 2 



i gnf ggrCO^ O) , „, 5 /,,-., 
5w d(vq) 2 



(£ a (ri,q;-l)(nn q ))a = -5— (C a (n, q; +1) (nn q ))n 
a z p a z p 



:L, (i + ^ £ )(.9t(0,0,0)~. 9t (0,^0)) 



5w d(vq) 2 

where gT(n,w,q) function is determined by Eqs. (IA4[) of Appendix lAl and 



|g I Up 3 q 2 i> 2 
3w 

<4<? 3 / flar(o,o,o) + 0gr(O,O,Oh 

5 a; \ dto 2 d{vq) 2 / 

|g|3,| fl*r(0,a,,0) ^ B 



5T (0,0,0) = 1 



d/y 



2/' 



From these results we immediately see that the correlation functions 7 a differ from unity only in the second order in 
v F \q\, i.e. 

la (q;P)^l + 0(T^ a a 2 P q 2 v 2 /oj 2 ). (69) 

From (|55aj) and (|68b|) we obtain that in the expression for the emissivity induced by the vector currents (f2U)) both 
scalar and vector components, (|67aP and (|67d|) , contribute at the order v F , 

4 ^ 4 v 4 

^(Xv,o(n,q)) ~ — eyap% T (0,w,0) , 
45 w 4 
2 ^2 4 

3(wXv,i («,'?)) ~ jT-lp eyap% T (0,w,0) . 

Working in the leading order in i>f we have to put "fv — > 1 in view of Eq. (|69|) . Thus, generalization of the corresponding 
result of Ref. [l8[ to arbitrary temperatures reduces in the leading order to the replacement 3g — » Qgx{0, u>, 0) with 

^(0,a, ) 0) = -_^= w i(l-2„(-)). (70) 
Finally, for the neutron PBF emissivity on the vector current we obtain (for one neutrino flavor) 



c 



nPBF 
vv,V 



AOn) 2 2 ±_ 4 2 2 _ -, /y, \ 



81 



This is precisely the result derived in Ref. [18|. The coefficient ^v F n distinguishes this result form that previously 
obtained in Ref. [l[ with the bare vertices, 

(0n) 4p G 2 A^ A ?^ 1 

_ 15tt3 l{ T h [ ) ~J ^TTi(e«w + 1) 2 ' 1 j 
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Expression (I7T1) deviates only slightly from the corresponding result obtained in Ref. [l6| with bare vertices. 

Now let us turn to neutrino emissivity induced by the axial- vector current. In the expansion v-p <C 1 the leading 
term contributing to the emissivity is of the order v^. Keeping only the leading terms we cast Eq. (|21|) as 

Ka « -9 2 A e 2 A pvl q 2 + % T (0, u , 0). (73) 

The correlation factors j a contribute at the sub- leading order ~ v™, therefore we neglected these terms in approximate 
expression (I75|) . We emphasize that the last two cross terms in the squared brackets in Eq. (|2ip cannot be eliminated. 
Finally, for the neutron PBF emissivity induced by the axial- vector current we obtain (for one neutrino flavor) 

^nPBF / i , 11 2\ 222 (On) fn A \ 

e uu,A - I 1 + 2i ~ 3 J 9 A e A v F.n tuu ' ■ (74) 

This again coincides with the result derived in 18]. The resulting emissivity is the sum of contributions (|7ip and (|74l) . 
We stress that Eqs. (|71[) and (|74l) are approximate expressions obtained in the leading order in vp. General result 
looks more cumbersome but it is easily recovered with the help of Eqs. (j67]). The latter equations are derived in the 
present paper at arbitrary temperature. 



V. SOLUTION OF LEGGETT EQUATIONS AND CORRELATION FUNCTIONS: AXIAL- VECTOR 

CURRENT 

In this section on example of the weak axial-vector current we will show how to correctly apply the Leggett 
formalism to calculate the current-current correlators. We operate with equations (|44|) and demonstrate explicitly 
how one should exploit the even and odd vertices (|38I39[) . 

The bare vertex generated by the axial- vector current (j8b[) gives rise only to the contribution to the t\ vertex. E.g., 
from (l2~5j) we have 

t"t = -QA— (vh - T) , t 1 =-g A (f Atl {n 1 q)l a -T Afl {n,q)l) , (75) 
a 

and the anomalous vertex is given by 

i i = -9A (TA,i(n, q) h - TA,o(n, q) f) . (76) 
In accordance with Eqs. (|38I39[) the even and odd vertices are 

*T" = ~9A — vlo, t^ = -g A (f AA {n,q) l - T A (n,q) I ) , 
ti U = 9A — 1, t 1 =-gA('r A!l {H,q)lo-T^ fi (n,q)l); 
T%i{n,q) = ±(f AA (n,q) ± T A ,i(-n,q)) , 

r A,o(™><?) = I ( r A,o(n, q) ±r A ,o (-ft,q)) ■ (77) 

Substituting these vertices in (|44[) and separating the parts proportional to aud I we arrive at two sets of equations 
for the scalar and vector vertices. The set for the t a and rj vertices is 

TA,o = rf ( - A t a , + A(r?', q)) n , r A ,o - (Y W '> 4,0 - <i </))„-, r^) , 

<o = V + r i ( - A K> *))*> ?U,o + <«(*', q)) 7r r£ + (4% k(u', q)),, r Afl ) , 

Taa = IT ( - A(n ', <?))., ta.o + K(n\ q)),, r+ + ((«(/?', <?) - 2A(n ', <?)))„- r Afi ) . (78) 

Here we explicitly take into account that the vertices do not depend on ft, provided the interaction constants and 
Tq contain only the zeroth Legendre harmonics. Since A and k are even functions of v q, these equations are simplified 

as 

ta,o = T i(- A ota,o + ( ~ X (n',q))^,T A fi - (fX(n,q))^, r A0 ) , 

T Xo = + r i ( K ( r? '' ?))*' r Xo> 

Ta,o = r i ( - ( 3 + (( K («''9) - 2A(tT', ?))>*, r7 )0 ) . (79) 
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The equation for t ao decouples from equations for t a and ta,i- Then the solution is 

r A ,o(n,q) =T+ fi (n,q)+T A0 (n,q) = t+ = ^ [l - Pf (/c(r?', g))^/] , T AQ = T Afi =Q. (80) 

Taking into account Eq. (| A30I) we have (/c(n , q))fi< — {L(n , q; — = {£ A (n , q\ — l))f?'j since (0{n. q; — 1))^ = 0. 
Thus, we recover Eq. (|54[) with Pa,o = — 1 ■ The anomalous vertex (I52[) vanishes for P A) o = — 1 in agreement with 

Eq. cm}- 

The second set of equations is for the vector vertices f A 1 and ta,i- 

TA,x{n,q) =rf (- (A T At i(n',q)) H > + ~ 2 % q} X(n' ,q)f A;1 (n' ,q))^, - (^\{n\q) f\i(n',q))^, 

~(fKn,q)T AA (n,q)) rr ^ , 
^Xi(™*,<7) = r i ((^ A(n',< ? )T j4a (n',g)) ? _,, + g) f ^(n g)), r + K(n , q) T Ajl {n , g))_.,J , 

r Ail (n,g) = Mff+ri- ((£ A(rr',«)^, 1 (rr',g)>, r , + (^/c(rT',g) f + 1 (rT', g))., 

+ ((K(n',q) - 2X(n',q)) f AA (n' , q))rt' 



The inspection of these equations reveals that the first and the third equations have trivial solutions f A 1 = ^ v and 
ta,i = . Presenting the solution for the second equation in the form f A1 = t+ n q we obtain 

t+ = Tf{{K{n\q)) n ,T+ + {K{n',q)U>^ 

= r i e ~t W\ W' I 1 - r i W> q)) n - 
In the last expression using first Eq. (IA30[) and then Eqs. (I61IIJ1[) we obtain 



— (K(n',q))fi' = (£ a (n',g;-l)) rT < 

\q\ \q\ 



(£ a (n',q;-l))fi> = (£ a (n',q;-l)(v'n g )}ri' 



(82) 



Thus, we recover expressions f|57l59[) . which we have derived above within the Larkin-Migdal formalism. 

Now we shortly dwell on the axial-vector current-current correlators, as they follow from the Leggett formalism. 
Substituting the vertices (|77p in Eqs. (|47p for the correlators and the corresponding polarization tensor terms we 
obtain 



xi 1 

xt 



Xi+ + xt 

2 e A ,2 / - 

a \ 

2 e A , J* / - 

a \ 

,2ttA 



jzn^K(n, q) rt i(n, q) + (/c(n,g) - 2A(n, g)) f A i(n, g) - ^-X(fi, q)f A . 1 {n,q) 
(vq) 'A 

K(n, q) r| (n, q) + {n{n, q) - 2A(n, g)) t a (n, q) - -^X(n, q) T A , (n, q) 



xt 



l' Il*(LO,q) + lo li Tl? (LU,q), 
-9 A — lk ( \^{n,q)r\ x(r?,g) + -y^r n(n,q) f A i(n,g) - X{n,q) T Ayl (n,q) 



i 2 ^ , • 
+ 9a~ 1 



K (n,q)TX fi {n,q) + y^r «(n, g) T Afi (n, q) - X(H, q)r Afi (n,q) 



(83) 



= z W>,g) + Ujn^,g). 

The elements of the polarization tensor IT^ derived here agree with those defined in Eq. ([T7]) with the correlators 
obtained within the Larkin-Migdal approach, (|d7|) . and the bare vertices (|8b|) 

We bring attention to the presence of the cross-terms IT^ and LT^, which were introduced in Ref. [T^]. These terms 
survive, even if we omit for simplicity all correlation effects and put = 0. In the latter case we get 



= Uf (u;,q) 



~9 A 



2 A 



a 2 \ (v q ) 



re(n, q) 



-9 A — — {K(n,q))n 
a q 



(84) 



Thus, we see that the cross-terms do not vanish, in contradiction with the claim of Ref. [22j]. The same cross-terms 
are obtained from Eqs. ©, ([TT]) . (152^1 and (p5"|) within Larkin-Migdal formalism. 
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VI. CRITICAL REMARKS TO REF. |2 



The results of Ref. [18| were criticized by Leinson in Ref. [22j, where the emissivity of the PBF processes in the 
neutron superfiuid with s-wave paring was recalculated in the framework of the Leggett approach. The Fermi liquid 
amplitude of the NN interaction considered in Ref . 1231 included the zero and first Legendre harmonics in the particle- 
hole channel. Comparison of the results of Ref. |22j with our results obtained in Ref. [l8| and here in Section IIVI 
reveals essential differences, even if one keeps only zcroth harmonics in the particle-hole interaction and also in the 
absence of such interactions. Therefore we are forced to do some relevant comments. 

(i) In Ref. [22| the spatial component of the weak vector current (the bare vertex oc v) does not contain the induced 
term oc q/\q\- It can be directly seen that such a term should appear in Eq. (31) of Ref. [22j . if we substitute T _ oc v 

in the third term of Eq. (31) and use the relation (|56p above. In the same way the anomalous vertex T gets the 
component oc q/\q \ . The presence of these induced components of the in-medium vertices is essential for establishing 
of the transversalit y of the polarization tensor and proving of the vector current conservation. Without these terms 
the Ansatz (39) in [22| is invalid. 

(ii) The temporal part of the axial- vector current (the bare vertex oc va) gets in Ref. [22| the renormalization factor 

l-fMn)-2A ( n)), * 1 + ^(0^,0) + 0(50 V * • (85) 

The result ([85]) is obviously incorrect, since the direction- independent interaction (zeroth harmonics) cannot renor- 
malize the vertex proportional to n, the loops in Fig. [3] do not depend on ft. 

From our calculations in Ref. [l8[ and here in Section[V]we see that only the average (fc(n))s ~ 0(v F q 2 /uj 2 ) enters 
the correlation factors in Eqs. (1501) and (|5Tj) . We can argue differently. We see that the first term on the right-hand 
side of Eq. (|85[) differs from unity even in the limit v-p — ¥ 0. Assume now uj 3> A. In this limit the correlation factor 
should take the form as in the case of a normal Fermi liquid with the ordinary Lindhard's function (taken now in the 
limit uj ^> \q\vp). It is well known that the Lindhard function is oc v F in this limit, see d 1, [26| and Eq. (|69j) above. 
This argument brought up in Ref. [HI was not perceived by the author of Ref. (22j ■ 

In Ref. |22| the author is surprised that in Ref. [l8[ we mentioned importance of correlation effects but dropped them 
in the final expression for the emissivity. We note that only in the expressions for the neutrino emissivities written 
up to the leading orders in vp (|71|73|74|) we can put j a — 1- This, obviously, does not mean that NN correlation 
factors can be always ignored. General expressions for vertices in Ref. [l8j contain correlation terms. Examples, where 
correlation terms give rise to important contributions, can be found in Refs. Also, in Ref. [Tij we mentioned 

a principal problem in applications of the BCS approximation to the description of pairing in nuclear systems. In 
the BCS approximation one uses the same interactions in particle-particle and particle-hole channels, whereas in the 
nuclear matter they can be significantly different. Since the Migdal theorem (valid for the electron-phonon interaction) 
does not hold in case of strongly interacting system, one should use the general Larkin-Migdal-Leggett formalism. 
Arguing, we mentioned that in particle-particle channel the Landau-Migdal parameter is necessary attractive /q < 
(to provide s-pairing), but in the particle-hole channel one has Jq > at some densities, as follows from estimates of 
these parameters, see [2tjj . 

The origin of the mistake that have led in Ref. [UJ to Eq. (|85|) is the wrong assignment of the bare vertex £± in 
Eq. (91) there. The physical vertex cannot be proportional to \v\ which in the coordinate space would correspond 
to the operator vr having bad analytical properties. Consequently, Eqs. (96) and (97) in Ref. |22| a re incorrect, 
as being written for scalars but not for vector objects. Following the original prescriptions of Ref7|27l HH, for the 
temporal part of the axial- vector current, v a, we have to put = v a and £_ = 0. Then the original Leggett 
equations (22) in Ref. [28[ yield the solution 



ipi = 0, -02 = ; , — , ^3=va. 



For the spatial component of the axial- vector current, a, we have £ + = and £_ = a and the solutions are 

ipi = $3 = , -02 = — j- -, zjz: — • 

1 - ( K (n ))*/ 

Thus, the main results of Ref. [22[ presented in Section VII and in Fig. 1 are false, since Eqs. (118,119) and (125) 
are derived with the incorrect solutions of the Leggett equations. 

(iii) The author of Ref. [22[ claims the vanishing of the mixed term of the polarization tensor n° l and IP , see Eq. 
(106) in [22j]. This conclusion is drawn at hand of incorrect Eqs. (98,99) in (22|. In Eq. (98) the left, bare vertex 
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must be taken as v and should stand under the integral, then the first term in the brackets gives non- vanishing result. 
In (99) the first term produces the finite value if one uses the results of the correctly formulated and solved Leggett 
equations. 

Further comments are in order: 

Ref. |22j | pays attention to the essential temperature dependence of the imaginary part of the retarded polarization 
function, which, as we have shown, is an artifact of the incorrect solution (|85|) . These corrections are small in the 
limit T <C A. Also they are small for ^ <C 1 for arbitrary T. In general, the imaginary part of the full retarded 
polarization function contains information not only on the processes of the one-nucleon origin (i.e. the PBF processes) 
but also on all multi-nucleon processes, e.g. on the two-nucleon NN- bremsstrahlung processes. It has been explained 
in Ref. [l|[ (see text before Eq. (8)). In the limit T <C 2A the phase space of different processes is well separated and 
contributions of the PBF and the NN- bremsstrahlung processes are easily decoupled. How to perform generalization 
to arbitrary temperatures was mentioned in Ref. (T^ |. also rough estimations of the dropped temperature dependent 
corrections are presented there, see discussion after (13). The full expressions for the emissivities of the PBF processes 
obtained in the given paper are obvious generalizations of the corresponding expressions of Ref. (l8l | . For «f < 1 the 
result of Ref. [18| proves to be correct in the leading order for all temperatures T < T c , provided for I(z) one uses 
general expression (|72[) rather than its low temperature limit. 

Even if we artificially omit all NN correlation effects, the result of Ref. [22| for the emissivity induced by the axial- 
vector current, Eq. (121), disagrees with the result of Ref. [l8j (see Eq. (34) in the latter work). The first difference 
is the factor (m*/m) 2 in the first term in the round brackets of (121), which Ref. |22j recovered in the form originally 
presented in Refs. Ref. [l|| argued that this term should be replaced by unity. The factor (m*/m) 2 appears, 

if one writes the temporal component of the axial- vector current as for the bare current ap/m = v m* /m, where v is 
the velocity of the nucleon quasi-particles. If the same arguments are applied to the spatial component of the vector 
current, there would be a problem with the conservation of the vector current. Indeed, the Ward identity between 
the Green's functions of the quasiparticles with energies e p = p 2 /2m* and the vertex is destroyed in this case (see 
a comment after Eq. (24) in [18|]). In general the central object of the study is the correlator of the fully dressed, 
in-medium currents and, therefore, the bare quantities cannot enter the final expressions. Appropriate excitations in 
Fermi liquids are quasiparticles obeying the Landau kinetic equation, where the group velocity v = (de/dp) \ PF = pf /rn* 
enters rather than the phase velocity pf /m. The proper Fermi liquid renormalization of the vertices was worked out 
in the seminal paper by Migdal [3^ . 

Ref. [22j stresses that the second correction term, 11/21, in Eq. (121) was obtained first in Ref. [![. Note that the 
correction found and used in previous works was twice as small, and it was corrected in Ref. [llj]. However the main 
difference between Eq. (121) in (22j) and Eq. (34) in [l8| is the absence in Eq. (121) of the cross terms of the temporal 
and spatial components of the axial- vector current, which falls out in Ref. [22j because of the incorrect solution of the 
Leggett equations. 

The factor j accom pan ying the Landau-Migdal parameter go in Eqs. (100,101) of Ref. [22[ with the origin in Eq. (90) 
was taken as in Ref. 1281 , where a stands for the operator of nucleon spin, whereas in the amplitude and the bare 
axial current in Ref. [22j |. see Eqs. (1) and (89) there, Oi enter as the Pauli matrices. 

In spite of the special commen t by Leggett that the relation (|43dj) is valid in the limit w|<f|,a; <C A only (see the 
comment before Eq. (17) in Ref. [28|V Ref. [22j uses this relation for oj > 2 A without any explanation. 



VII. CONCLUSION 



In Ref. [25| Larkin and Migdal extended the Fermi liquid approach onto Fermi systems with pairing. The equations 
for the full normal and anomalous vertices (the Larkin-Migdal equations) have been derived. They considered a 
particular case of s-wave paring at zero temperature, aiming at applications to atomic nuclei {33T - I351 ] . In Refs. [27H281] 
Leggett generalized this approach to the Fermi liquid at non-zero temperature and applied it to study the low- 
frequency, low momenta collective excitations. In difference with Larkin and Migdal, Leggett formulated equations in 
a matrix form for symmetric and antisymmetric vertices (the Leggett equations). Explicit equations were formulated 
for vertices with the symmetry (1, cr 3 , p^, p^cra) . 

The present analysis provides necessary, although self-evident, extensions of the Larkin- Migdal formalism (25j to 
arbitrary temperature and of the Leggett formalism [271 . |28| to arbitrary frequencies and momenta and to the vector 
and axial-vector weak current symmetry Efficiency of the Larkin-Migdal and the Leggett approaches generalized in 
such a way is demonstrated on example of the calculation of the neutrino emissivity in the reactions of the nucleon 
Cooper pair breaking and formation. To be specific we considered s-state neutron pairing and included only zero 
harmonics in the Fermi liquid interaction both in the particle-particle and in the particle-hole channels. Compared 
to our previous work [l8j . where explicit expressions were found in the low temperature limit T <C 2 A, here we 
performed generalizations to arbitrary temperatures. Since recently there were published works, where presence of 
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the formal difference in the Larkin-Migdal and Leggett approaches have led to misleading conclusions, we carefully 
analyzed both approaches demonstrating explicitly that they indeed lead to the very same results. 

First, from the diagrammatic equation for the current-current correlator and vertices we reproduced the Larkin- 
Migdal and the Leggett equations. Within the Matsubara techniques for T^Owe showed the correspondence between 
the vertices and the loop functions introduced in both approaches. It turned out possible to cast all necessary loop 
functions in terms of the one master function qt- In passing we proved the validity of a relation between the loop 
functions at arbitrary frequencies and momenta, which has been previously proven by Leggett in the low frequency- 
momentum region. 

Then we solved first the Larkin-Migdal and then the Leggett equations for the renormalized vertices of the neutral 
weak currents at arbitrary temperature and found the current-current correlator for the neutrino-antineutrino pair 
in the neutron superfluids, imaginary part of which describes the neutrino-antineutrino pair production in reactions 
with the nucleon Cooper pair breaking and formation. Also we proved the exact conservation of the vector current. 

The results of recent publication [22J based on the Leggett approach are found to be invalid because of falsely 
interpretation of the Leggett equations and their wrong solution. 
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Appendix A: The loop functions 



At zero temperature the loop functions (|3~4"|) were calculated in Ref. [3(j using the Feynman method for the integral 
of the Green's function products 



L(n,q;P) = a 2 p 

M(n,q) = - 

N(n, q) — a 2 p 

0(n,q;P 

The universal function g is given by 

g = 2 / d$ 



qy n , \ \ .9(2) n , m 
_ w (i-,(,))-_(i + P) 

2 u + qv 
a P oa 9{z), 



2A 

uj 2 — (qv) 2 



9{z), 



a 2 p 



4 A 2 
uj + qv lu — qv 



4A 



4A 



(Al) 



F,F_ 



a 2 p 



2 A 2 d$ 



[e\-E\] {e 2 _-E 2 _}> 

where J d$o = / d$T for T = 0, see Eq. (|2"T)) ; e± = e ± \w and E± = E p ± q f 2 . Calculation of the integral yields 



9(z) = 



arcsinh\/ z 2 — 1 iird(z 2 — l) 



zx/I 2- ! 2z N /i 2_ T ' 

uj 2 - (qv) 2 

> 1 , i) = %ti. 



(A2) 



4 A 2 

At finite temperatures the Feynman method does not work [36| and the Matsubara techniques can be used instead: 

+00 

rp 

" (A3) 



d$ r ^- = 2A 2 / de p ^ —— 

a P J n=-oo ^ len 



-E 2 + ) ^e n ) 2 -E 2 _Y 



Here e n = (2 n + 1) tt T and u) m — 2 m tt T . Thus, we obtain 



de r , 



(E+-E-) (n(E-) - n(E+)) (E+ + E-) (1 - <g_) - n(E+)) 
E+E_ (iuj m ) 2 - (E+ - E_) 2 E + E_ {iLu rn ) 2 - (E+ + E_) 2 



(A4) 
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with the fermion occupation function n(x) = l/(exp(x/T) + 1) . After the replacement iuj m — > uj + — uj + i we obtain 
the analytical continuation to the retarded <?t function in the uj complex plain. Other Matsubara sums of products 
of the normal and anomalous Green's functions can be obtained with the help of the following general relation 



where 



S ab (e p , n, iu) m q) = T ^ Jf- 
= (n(E_)-n(E + )) ( 



va 



X. 



[(»£„ + iu m ) 2 - E\] [(ie n y - El 

yayb 



YlY b _ 



iui m — E + + E- iu> m + E + — E- 



+ (l-n(E.)-n(E+)) (- 



Y + Y1 



Y?Y*_ 



iu m + E + + £!_ iuj m — E+ — E- 



(A5) 



B 



ie n + iui m + e+ , B n _ = ie n + e_ 



B n ,+ = ~ie n - iu m + e+ , B n _ = -ie n + e_ , 

Y± = {A,A±,A ± y2E ± , Y± = {A,A ± , A ± }/2E ± , 

A± = e± + E± , A± = e±-E±, e± = ep±$/2 ■ 

Convolutions of the Green's functions can be expressed through the corresponding elements of the matrix S ab inte- 
grated of de p , e.g., 

{F + -G h _){n,uj,q) = a 2 p [ de p S 13 {e p , n, uj, q ) . 



The matrix (| A5|) has the following properties 



and since Y^ ] = Y*± 



S ab (e p , n, uj, q) = S ba (-e p , n, -u, -q) , 



S (e p ,n,uj,q) = S (e p ,n,-ui,q) 



(A6) 



(A7) 



Relation (|A6|) allows to interchange the order of Green's function in the convolutions with the simultaneous change 
of uj — > — uj and q — > —q, e.g., 

(F+.G h _)(n,UJ,q) = (G\-F-)(n, -uj, —q) , 
(G+-G h _) (n, LU,q) = (G h + -G-)(n,-u,-q), 

(F+-G-)(n,uj,q) = (G+-F_)(n, -u, -q) . (A8) 

From (|A7|) we find 

(G h + -G h _){n,u>,q) = (G + -G-)(n, -uj,q) . (A9) 
Now we are able to consider the Green's function products entering the matrix (J3BJ). We find, see fill. [l7j|. 

(G + -G-)(n,u>,q) = a 2 p / de p 



n(E-) - n(E+) (e+ £_ + E+ e_) w + (e+ e_ + £+ E_){E+ - £_) 
2E+E- w 2 - (£+ - £_) 2 

1 - n(E-) - n(E + ) (e+ £L - e_) w - (e+ e_ - £L) (J5+ + EL) 



2E+E. 



^-{E + + E_y 



, (A10) 



+00 



(G+-F_)(n,o;,g) - -« 2 pA / de 



n(E_) - n(E+) ojE + + e+ (£L - E_ 
2E+E- uj 2 -(EL -EL ) 2 

OO 

1 - n(£_) - n(E + ) -cj£; + - 6+ (£+ 



2E+E. 



uj 2 -{E++E^f 



(All) 
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+oo 



(F + -G-)(ri,u>,q) = -~a 2 pA / de. 



n{EJ) - n(E + ) w£L + e_ (E + - E_) 



+ 



2E+E- lo 2 -(E + -E-) 2 

OO 

1 - n(E-) - n(E + ) w - e_ (E+ + E_ 



2E + E_ 



lo 2 -(E + + E_) 



(A12) 



Making the replacement e p — > — e p in the integral (| Al 1|) . which induces the changes e± — > — e T and E± — > £? T , we 
obtain exactly the same integral as in Eq. (|A12I) but with opposite sign. Thus, we prove that 



(G + -F_)(n,u,q) = -(F+-G-)(n,u,q). 



(A13) 



For the (G+-G^) loop we have 



+00 



(G+-Gi)(ft,u,$) = -a 2 p / de 



+ 



n{EJ) - n{E+) (e_ E + - e+ £L)w + (e+ e_ - E_ E+) (£L - E+) 
2E+E- uj 2 - (E+ - E-) 2 

1 - n{E-) - n(E+) -(e_ E+ + e+ E_)lj + (e_ e+ + E+ £L) [E+ + E-) 



2E+E- 



uj 2 -{E+ + E_) 2 



(A14) 



It is easy to verify that this integral does not change under the replacement q — > —q, since then e± — > e T and E± — > E T , 
and therefore (G + ■ G h i)(n,uj,q) = (G + ■ G h _){n 1 uj 1 ~q) . One can see that the replacement e p — > — e p in the integral 
is equivalent to the replacement to — > —to. Hence we prove the relation (G+ ■ G^)(n, — to,q) = (G+ • G ft i)(n, to,q) . 
Combining the last two equations with (IA8|) we prove 

(G+-G'l){n,uj,q) = (G+-G h _)(n,-uj,-q) = (G h + • G_)(n, to, q ). (A15) 

Now we state the results for {F + -G l t) and F-) products 



+00 



(F + -G h _)(n,uj,q) = ~a 2 pA / de p 



n{E-) - n{E+) -E- uj + e_ (E+ - £_) 



2E+E- u 2 -{E+-E-) 2 

— OO 

1 - n(E-) - n(E+) -E- to - e_ (E + + EJ)~ 



2E+E. 



u 2 + (E+ + E_y 



(A16) 



+OO 



(G5.-F_)(fJ,w,g) = -a 2 pA / de 



ra(£L) - n(E + ) -E + lu + e+ (E + - E-) 
2E+E- u: 2 - (E+ - E-) 2 

— OO 

l-n(-EL) - n(E+) E+io - e+ (E+ + EJ)- 



2E+E- oj 2 + (E + +E-) 2 

Making use of the replacement e p — > — e p as in Eqs. (|A11IA12|) we are able to prove that 

{F + -G h _){n,u,q) = -(G h + -F_)(n,to,q). 



Eqs. (|A13I) . (|A15I) and (|A18|) are used in the paper body, see Eq. (|28]l, 

Now we turn to the derivation of relations (|43|) . Consider first (|43aj) . From (|A12j) and (|A17|) we have 



+00 



(G+-F-) - (F + -G-) = -a 2 P A / de 



n(£L) -n(E+) -{E + +E_)lu + (e+ - e+) (E+ - E-) 



2E+E_ 



+ 



uj 2 — {E + - E-) 2 

— OO 

1 - n(E-) -n{E+) (E+ - £L) to - (e+ - e_) 



2E+E- 



to 2 + {E++E-) 2 



(A17) 



(A18) 



The terms with to in the numerator vanish exactly since they are antisymmetric with respect to the replacement 
e p — > — e p . In the remaining terms using e+ — e_ = vq and comparing with (|A4|) and (|A3j) we find 



(G^-F_) - (F+-G-) = -a 2 p^Lg T (n,uj,q-) 



v q 



(F+-F_). 
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Similarly, using the symmetry properties of the integrand under the change e p — > —e p we verify that in the sum 
(G+ • F-) + (F + ■ G ) only the term with oj will survive in the numerator 



+00 



(G^-F_) + (F+-G_) = -a 2 pA / de p 



n(E-)-n{E+) -(E + -EJ)oj 



2E+E- uj 2 -(E + -E-) 2 

— OO 

l-n(E-)-n(E + ) (E + +E_)oj 



2E+E- uj 2 + (E + +E_) 2 
= a P ^ 9T(n,uj,q) = — (F+-F-) . 

Thus, we recovered Eq. (|43b|) . 

Before we consider relations (I43c|) and (|43d|) let us simplify the numerators in (IA10[) and (| A14|) . We use 



and obtain 



e± = e p ±\(vq), + E 2 _ ) = e* + \ (vq) 2 + A 2 , \ (E 2 + - E 2 _ ) = e v (vq) , 



(e+ E_ ± E+ e_) oj - (e+ e_ ± E + EJ) (E- T £+) 

= (£L T £+)(a 2 + - (25- ± ^ + ) 2 )^SP 



(A19) 



(A20) 



From Eq. (| A9|) follows that (G\ ■ G 1 !) — (G+ • G-) is an odd function of oj and contains only the terms from (|A20p 
linear in oj. Oppositely, (G+ ■ G_) + (G+ • G_) is an even function of oj and contains the terms from (|A20[) independent 
of oj. Hence we can write 



+00 



n(£L) - n(£ + ) _ ) 2 - (#+ + ^-) 2 



2E+E_ 



1 -n(E-) ~n(E + ) 
2£ +J E_ 



w 2 -(£ + -£_) 2 
vg) 2 ) 2 " 



w 2 - (J5+ + £_) 2 



(A21) 



and 



• G h _ + G+ ■ G_) + (F+ • F_) = a 2 p \ I de 



- E+) - 



2E+E_ 



1 - n(£L) - n(.E7 + ) 



w 2 - (£ + - 
(z/g) 2 ) 2 - 



u; 2 - (£7+ + 



Comparing Eqs. (|A21I) and (|A22|) we conclude that 



(vq) 



\ (G\ ■ G h _ + G + ■ G- ) + (F+ ■ F_ 



(A22) 



(A23) 



thus, recovering Eq. (|43d[) . 

We turn now to the derivation of the relation (|43c[) . First, we note that in (|A14j) the terms with oj in the numerator 
do not contribute because of their symmetry properties with respect to the change of the sign of e p . In the other 
terms in the numerator we can use 

e+ e_ ± E+E- = \(E+ ± £L) 2 - \ (vq) 2 - A 2 . 

Then after the separation of the pole part we obtain 



(G+-G h _) = a 2 p 



oj 2 - (vg) 2 - 2 A 2 1 
2A 2 2 



g T (n,oj,q) 



-\--yo 



+ a z p- 
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QJ 2 -(vq) 2 -2A 2 1 
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+00 



gr(n ,oj,q) + a 2 p I de 



1-2 n(E p ) 

2 En 
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The last divergent term must be cut of at some scale £ ~ p and is nothing else but the quantity Aq in Eq. ([5]), which 
enters the gap equation @. The relation (|43c|) follows now immediately. 

Now we are in the position to generalize expressions for the loop functions (|34j) for T^O: 

M(n,u),q) = -a p g T (n,uj,q), 



_ 2 lu 2 - (qv) 2 , _ 
N(n,ui,q) = a p — -g g T (n,uj,q), 



0(n ,ui,q; P) = a 2 p 



4 A 2 

w + (ft/ cj — qv 



P 



4A 4A 

with gr given by Eq. ()A4|) . The L function is less straightforward. We find 



g T (n,uj,q] 
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(1 - n(EJ) - n(E + )) 
oj 2 -(E++E^) 2 

Then with the help of the relation 



de r . 



2E+E. 



(n(E_)-n(E + )) 



.uj 2 -(E + -E-) 2 

— OO 

(E + + E_)((vq) 2 -(E_-E + f) 



(E + - E_)((vq) 2 - (E_ + E + ) 2 



(vq) 2 - (£_ ± 



4(vq) 2 A 2 



(E- TE+) 2 - (vq) 2 



we present 



(vq) 2 - (£?_ ± £+) 2 4 (tTg ) 2 A 2 



oj 2 - (E+ T E-) 2 uj 2 - (vq) 2 \uj 2 - (E+ q= £-) 2 («7) 2 - (£L T £+) 2 

Finally after comparison with Eq. (|A4|) we derive 

L(n , ui,q;P) = a 2 p 



Vq ^^ (gr(n, (vq),q) -g T (n,ui,q)) - 1 ~|~ g T (n,ui,q) 
u> — v q x '2 



(A25) 



(A26) 



(A27) 



(A28) 



For T = 0, gr(n , (vq), q ) = 1 and the old result (|A1|) is recovered. 

The relation between the function n used in the Leggett equation (|44"1) and the functions used in the Larkin-Migdal 
equations (j3~4")l is derived with the help of the substitution of Eq. (|A27|) into Eq. (IA22[) and by making use of the fact 
that the function gx is even function of q. Thus we obtain 



K(n,uj,q) = a 2 p 



ui 2 — (v q) 2 



(g T (n, (vq),q) - g T (n,u),q)) 



v q 



lu + v q 



L(n,u,q;-1) = \ (L(n ,ui, q; -1) + L(n, -uj, q; -1)) 



(A29) 



From (|A28[) and (|A4[) we verify that L(n , — uj, — q; P) — L(n,ui,q;P) and L(—n,ui,—q;P) — L(n,ui,q;P), and 
therefore 



(K(n,uj,q))n = \ ((L(n , uj, q ; + (L(n , -uj, q ; -l))n) 
= \({L(n,u>,q;-l)) R + (L(-n , uj, q ; 
= {L(n,ui,q;-l))n- 



(A30) 
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Appendix B: Current conservation 



Here we prove the transversality of the polarization tensor (jl3|) , IP" q v — which guarantees the conservation of 
the vector current in a superfluid Fermi liquid. 

As we have shown in Appendix \^ the loop functions M, N, O and L are expressed through the only g? function 
defined for arbitrary temperature. From definitions ([55]). (|60|) of functions C a and C a and expressions for the loop 
functions (|A24I) and (|A28j) we deduce that 

lo (C(n,q;P afi ))n = (£(n, q; P a ,i) {vq))a ■ (Bl) 

From Eq. (fl"3|) we obtain 

nv*(w,?)g,i = — w(xy,o {n,q))n- (qxv,i(n,q))n , (B2) 

ny(w ;( f)(7 M = ^ (wxy,o(n - ?Xf,i("'.9)))b . ( B3 ) 

with ey = 1 . To demonstrate transversality of the polarization tensor for the vector current we need to prove 
vanishing of these components. 

Applying the averages given in Eqs. (|67al67cp . Eq. (|B1[) and using that P aj o = —Pa,i we obtain 

(wXa,o(n, q) -qXa,i(n, q))n 

= W7 a (g; P„ )0 ) (C(n, q; P a ,o))n T a,o ~ Ja(q; -P a ,i) {£(ft, q; P a .i) (n • q)) n 
= <o7a(g; Pa.o) (w (C{n, q; P a , )) - (£(n, q; P a ,i) (v ■ <f ))) = . 

Remarkably, this relation holds both for the vector current and for the axial- vector currents, i.e. n° M (w, q ) q^ = 0. 
Now we turn to the second necessary relation (|B3[) . Calculating the product qxv,i we use Eq. (JS5J). Consider first 

{{uj 1 - {vq'Y)g T {n,q))-, 
+ lv {q- +1) r v q- +1) (£y(< <?; -l)(gn')) ff - C v (n, q; -1) (gn) (£ v (n, q; t£ x ■ 

Here we used Eqs. (|A24I) . Making use of this relation and Eqs. (IBTj) and (|B1[) we calculate 

((u<f) (wxv )0 - qxv,i)) n 

= (fo«?) [ujv(q;+~L)£v{n,q;+l) - (vq)j v (q;+l) £ v (n,q;~l) - q8xv,i{n,q)]) n Ty fl 
= {({vq)C{n,q-+l)) n uj - ({v qf C(n,q;-1)) H ) r£ 

((vq)(uj + qv)g T (n,q))f i 2 

~ a P-n—2 — T^xT\ h \\ (fog) 9T{n,q))n- (B4) 

((W - {yqy)g T {n,q))ti 

Interestingly, there is no dependence on Tq in the last expression. Finally, using explicit expression for the loop 
functions (|A24IA28|) we arrive at the identity 

(vi (wxv,o(n,gO - qxv,i(ft, q))) n = e v a p{v .(v q))n , (B5) 
which means that Tly (w, q) q^ = 0. 
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